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Abstract 

Nonlinear Schrodinger equation with the Schwarzian initial data is 
important in nonlinear optics, Bose condensation and in the theory of 
strongly correlated electrons. The asymptotic solutions in the region 
x/t = 0(1), t — > oo, can be represented as a double series in t^ 1 and 
Int. Our current purpose is the description of the asymptotics of the 
coefficients of the series. 
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1 Introduction 

A coupled nonlinear dispersive partial differential equation in (1 + 1) dimension 
for the functions g+ and 

-id t g+ = \ d l9+ + ^glg-, 

idtg- = ^d 2 x g-+4g 2 _g+, (1) 

called the separated Nonlinear Schrodinger equation (sNLS), contains the con- 
ventional NLS equation in both the focusing and defocusing forms as g + = 
or g + = — g_, respectively. For certain physical applications, e.g. in nonlin- 
ear optics, Bose condensation, theory of strongly correlated electrons, see JjJ 
- ||, the detailed information on the long time asymptotics of solutions with 
initial conditions rapidly decaying as x — > ±oo is quite useful for qualitative 
explanation of the experimental phenomena. 

Our interest to the long time asymptotics for the sNLS equation is inspired 
by its application to the Hubbard model for one-dimensional gas of strongly 
correlated electrons. The model explains a remarkable effect of charge and spin 
separation, discovered experimentally by C. Kim, Z.-X.M. Shen, N. Motoyama, 
H. Eisaki, S. Ushida, T. Tohyama and S. Maekawa p9|. Theoretical justification 
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of the charge and spin separation include the study of temperature dependent 
correlation functions in the Hubbard model. In the papers p]]-||, it was proven 
that time and temperature dependent correlations in Hubbard model can be 
described by the sNLS equation 

For the systems completely integrable in the sense of the Lax representa- 
tion |[o| the necessary asymptotic information can be extracted from the 
Riemann-Hilbert problem analysis Often, the fact of integrability implies 
the existence of a long time expansion of the generic solution in a formal series, 
the successive terms of which satisfy some recurrence relation, and the leading 
order coefficients can be expressed in terms of the spectral data for the associ- 
ated linear system. For equation ([!]), the Lax pair was discovered in while 
the formulation of the Riemann-Hilbert problem can be found in |^| . As t — ► oo 
for x/t bounded, system (|l|) admits the formal solution given by 

n=l k=0 

oo 2n A,\k 

9- = e-**r-Ci-*0 * « (, + E E ^**) > (2) 



where the quantities v, uq, vq, u n u and v n k are some functions of Ao = —x/2t. 

For the NLS equation (<? + = ±<?_), the asymptotic expansion was suggested 
by M. Ablowitz and H. Segur ||. For the defocusing NLS (5+ = —gJ), the 
existence of the asymptotic series (|^) is proven by P. Deift and X. Zhou || 
using the Riemann-Hilbert problem analysis, and there is no principal obstacle 
to extend their approach for the case of the separated NLS equation. Thus 
we refer to (||) as the Ablowitz-Segur-Deift-Zhou expansion. Expressions for 
the leading coefficients for the asymptotic expansion of the conventional NLS 
equation in terms of the spectral data were found by S. Manakov, V. Zakharov, 
H. Segur and M. Ablowitz, see O] [[[g). The general sNLS case was studied 
by A. Its, A. Izergin, V. Korepin and G. Varzugin [[17 |, who have expressed the 
leading order coefficients Uo, vo and v = —uqVo in (|2) in terms of the spectral 
data. 

The generic solution of the focusing NLS equation contains solitons and 
radiation. The interaction of the single soliton with the radiation was described 
by Segur fl8|| . It can be shown that, for the generic Schwarzian initial data 
and generic bounded ratio x/t, \c — wz\ < M, c = const, where M is small 
enough, the expansion (^) for sNLS equation (Q) represents the contribution of 
the continuous spectrum when possible solitons go away from the moving frame. 
Thus the soliton contribution in the considered asymptotic area is exponentially 
small for large t and can be neglected in compare with contributions of the 
power-log terms. 

Below, we pay our special attention to the coefficients u n> 2m ^n,2m which 
determine the leading behavior of g + and g_ as t — > 00 in the order t~ n . For 



2 



these coefficients as well as for u n ,2n-i> v n ,2n-i, we find simple exact formulae 

i n (v') 2n {-i) n {v') 2n 

Un,2n=U — , V n>2n = V — { , (3) 

and ( pp[ ) below. We describe coefficients at other powers of In t using the gener- 
ating functions which can be reduced to a system of polynomials satisfying the 
recursion relations, see (p4|), (|23| ). As a by-product, we modify the Ablowitz- 
Segur-Dcift-Zhou expansion (0)7 



9+ = exp 



g- = exp 



. x 



1 1 . A .(^) 2 ln 2 4ti ^ ' (ln4t) fc . 



. 1 . , V 111 11 v — ■( v — ■( 1 111 1 

/ — - (- + %u) In 4* + j ^T-' ; " 



2i - 

n=0 k=0 
n=0 fc=0 



2 Recurrence relations and generating functions 

Substituting (||) into (Q), and equating coefficients of t^ 1 , we find 

v = -u v . (5) 

In the order t~ n , n > 2, equating coefficients of In 3 4t, < j < In, we obtain 
the recursion 

iv" (v') 2 

— + l)U n ,j+l + iniln.j = VU n j — U n -X,j-X Z — u n-l.j-2 - 

o o 

iv' 1 n 

— 4" u n-l,3-l + o U n-l,j + E E u i,««k,/?Um,<y, (6) 

t,fc, 777 = ce=0,...,2i 
Z + fc + m = ji /3 = 0,...,2fc 

1=0 Jm 

a + /3 + ~,=j 

iv" iv') 2 
o o 



IV . 1 ., x - \ - 



Ul, a Vk,l3V m ,-y, (7) 



777=0 a=0,...,2i 

l+k+m=n ^ = 2fc 

7 = 0,. ..,2777 



where the prime means differentiation with respect to Ao = —x/(2t). 

Master generating functions F(z, C), G(z, C) for the coefficients u n ^, v nt k are 
defined by the formal series 

F(z,C)=5>„,^"C fc , G(^,C)-^^^"C fc , (8) 
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where the coefficients u nt k, «n,fc vanish for n < 0, k < and k > 2n. It 
is straightforward to check that the master generating functions satisfy the 
nonstationary separated Nonlinear Schrodinger equation in (1 + 2) dimensions, 

-iF< + IZ F Z = {y- - { ^fzt 2 )F - *£zCF> + \zF" + F 2 G, 

iG c - izG z = (v+ - ^< 2 )G + %z<@ + \zG" + FG 2 . (9) 

8 8 4 8 

We also consider the sectional generating functions fj(z), gj(z), j > 0, 

oo oo 

fj( z ) = X! u n.27i-jZ n ) 9]( z ) = X! v n,2n-jZ n - (10) 
n=0 n—0 

Note, fj(z) = gj{z) = for j < because u n ^ — v n ^ = for k > 2n. The 
master generating functions F, G and the sectional generating functions fj, gj 
are related by the equations 

oo oo 

Using ( pi] ) in (^|) and equating coefficients of we obtain the differential 
system for the sectional generating functions fj (z) , gj (z) , 



f w " f f w ' , i f 4- ft 

fc + i + m =j 

2izd z gj-i - i(j - - = (12) 

if" (^) 2 w 1 , l „ J-y . 

= v 9j + z —9i-i ~ z ~ ^— 9j + z —9 3 -i + z g9j-2 + 2^ Jk9i9m- 

k,l,m=0 
k + l + m=j 

Thus, the generating functions fo(z), go(z) for u n .2n, v n ^ n solve the system 

izd z f = vf - z^-fo + f 2 g a , -izd z g = vg Q - z^-g + fag 2 . (13) 

The system implies that the product fo(z)go(z) = const. Since /o(0) = Uq and 
ffo(0) = voi we obtain the identity 

f g (z) = -v. (14) 

Using (O) in jl^), we easily find 

f (z)=uoe*— z =u }2-^ 

n=0 



-Z 
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-i^z ^ {-i) n {v') 2n „ 
9o(z)=v e « * = v J2 y g» ' *' ( 5) 

n=0 

which yield the explicit expressions (||) for the coefficients w nj 2n, Un,2w 

Generating functions fi(z), gi(z) for u„.2n-i, v n ^ n -i, satisfy the differential 
system 

II/" (l/') 2 ^ , 2 

-2tzd z f + izd z fi = vfi - z —fo - z ~ ^— h - z ~fo + 2 .hfo9a + f gi, 

iv" iv') 2 iv' 

2izd z g - izd z gi = ug 1 + z—go ~ z ~^~9i + z —9o + /i5o + 2 /o.9o.9i- (16) 

We will show that the differential system ( [l6| ) for f\(z) and g\(z) is solvable in 
terms of elementary functions. First, let us introduce the auxiliary functions 

fo(z) g a {z) 
These functions satisfy the non-homogeneous system of linear ODEs 

^ iv i -> i(y') 2 v" v' g' Q 

so that 

(v 2 )" 

»*(g l +Pi) = -V-' ( 18 ) 

oV 

where we have used the relations ( |l4|) an d (|l5| ). Taking into account the initial 
condition pi(0) = ?i(0) = 0, the system (|17|) is easily integrated and we arrive 
at the following expressions for the generating functions /i, gi, 

fi(z) =Pi(z)f (z), f Q (z)=u e l « z , 

, , , ivv" v" v'u' (v'fv" - 

»i(z) = ( -Iz - i- — z , 

F y ' V 4 8 4u ; 32 

Si(*0 = Qi{z)9o («)) ffo(*0=Voe 1 8 z , 

01(2) = ( ?-)z + i^- ■ z 2 . (19) 

yu ; V 4 8 4w ; 32 y ' 

Finally, expanding the generating functions in series of z, we find the coefficients 

«n,2n-l and v n ,2n-l, 

, ivv" v" v'u'n\ Aw" v" v'v'n^ 

^ = U ^-— ~ ¥ - ^ = - T - 4^ (20) 
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lt»,2n-i = -2n g „; ' s, y" + 7 + ivy" + -) , n > 2, 

8™(n — 2)! V n — 1 2 zto / 

Wn.an-X = -2^o V V ' , + r ^ ~ ™» + ~ ~) ) , " > 2 - 

8 n (n — 2)! V n — 1 2 vq / 

Generating functions fj(z), gj(z) for u n ^ n -j, Vn,2n-j, j > 2, satisfy the 
differential system (^2j). Similarly to the case j = 1 above, let us introduce the 
auxiliary functions Pj and qj , 

Pi = 4, 3i = ^- (21) 
Jo So 



In the terms of these functions, the system fll2|) reads, 



d z p 3 = — (Pj +qj) + dj, d z qj = — + q 3 ) + bj , (22) 



where 



aj = 2d zPj ^i + (— — — )Pi-i - 

y' (Pj-ifo)' .(Pi-2/0)" . ^ 

6j = 20,<&_i + ( — - —^-)Qj-i - 

y' (qj-igoY , , (gj-2go) ;/ w 
--: M — 3 > Pkqiq m - [to) 

4 3o 8,70 « fc ,f^ =0 

With the initial condition ^(0) = qj(0) = 0, the system is easily integrated and 
uniquely determines the functions Pj(z), qj(z), 

Pj (z)= [' aj (Qdc+iu f ^ / de( aj -(o+6i(0) 5 

Jo Jo s Jo 

?iW = f 6j(0 rfC -™\ ZA T f C df(ai(0 + 6i(0)- (24) 
Jo Jo s Jo 

These equations with expressions ( p^ ) together establish the recursion relation 
for the functions Pj(z), qj(z). In terms of Pj(z) and qj{z), expansion (Q) reads 

+ ^ ln J 4t 

3=0 



9 _ =e -^-(|-w)h«-^^^M^E), (25 ) 



j=o 



G 



Let cij(z) and bj(z) be polynomials of degree M with the zero z = of 
multiplicity m, 

M M 

k—rn k—m 

Then the functions Pj(z) and qj(z) ( |2~i| ) are polynomials of degree M + 1 with 
a zero at z = of multiplicity m + 1 , 



fc— m+1 
M+l . . 

*( z )= S l(hk-i-j(a 3 ,k-i+b^ k -i))z k . (26) 

fc— m+l 

On the other hand, cij(z) and 6j(z) are described in p^ ) as the actions of the 
differential operators applied to the functions py, qy with j' < j. Because 
Po(z) = Qo(z) = 1 and pi (z) , Qi(z) are polynomials of the second degree and 
a single zero at z = 0, cf. (|19|), it easy to check that 02(2) and &2(z) are non- 
homogeneous polynomials of the third degree such that 



(V) 4 (VM 2 (v'Yiv") 2 
«2,3 = - ^ y j (2-iv), 6 2 , 3 = - 1 j 2l l J (2 + (27) 



1/ f Mg «Uq , IVV V 



_ - "U -"U L 

fl 2,o — : 7, : 7: — j 02,0 



<± o 4«o 8uq ' 4 8 4«o 8t>o 

Thus P2(z) and 92 (z) are polynomials of the fourth degree with a single zero at 
z = 0. Some of their coefficients are 

p^^-MVT!, (28) 



4 v 4 8u 4wg 

-<.-*,)£ -*)-"«>" 



2 ' 



4 v 4 8u 4u, 
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Lemma 2.1 Functions Pj(z), Qj(z), j £ Z+, defined in (21) are polynomials of 
degree not exceeding 2j . The zero of pj(z) and Qj{z) at z = has multiplicity 
no less that mj — [^-1 • 

Proof. The assertion holds true for j = 0, 1, 2. Let it be correct for Vj < j'. 
Then aj'(z) and bji(z) are defined as the sum of polynomials. The maximal de- 
grees of such polynomials are deg((py_i/ )7/o) = 2j'-l, deg((qy_i0o)7ffo) = 
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2f - 1, and 

J 1 - 1 . .-, J"- 1 



z / \2 



Thus degaj'(z) = degbj'(z) < 2j' — 1, and degpj>(z) — degqj/(z) < 2j'. 

Multiplicity of the zero at z — of ay (z) and by (z) is no less than the min- 
imal multiplicity of the summed polynomials in ( |23|) , but the minor coefficients 
of the polynomials 2d z , py^\ and —(J — l)pj<_i/,z, as well as of 2d z qy-\ and 
—(J — 1)<7j'_i/z may cancel each other. Let j' = 2k be even. Then 

my = min{m 3 -'_i; my-i + 1; min (m a + mp + m 7 ) } = ^— — [— ^ — ] . 

Q + /3 + -/=3' 



Let j' = 2k - 1 be odd. Then 2my-\ - (j' - 1) = 0, and 
my = min{TOj'_i + 1; my _2 + 1; min (m a + + m 7 ) } 

a,0,"f=O,...,j f — 1 
Q + /3 + 7=3' 



/ + 1 



□ 

We arrive at the following expansions of the polynomials Pj(z), Qj(z), 



The use of (gg) in (|25|) yields (^). On the other hand, using (|29|) in (|2J), we 
obtain the expressions for the generating functions fj(z), gj(z), 

t ( \ n l k {v'f k 

f j {z)=u Q ^ z ^ p j|W _ fc gfcfc! , 

n=[i±i] fc=max{0;n-2j} 

»w = «o 2^ 2 2^ — gT^ — ' ( 3 °) 

n=[i±i] /c=max{0;n-2j} 

and the coefficients of our interest, 

i k {v') 2k 



k— max{0;n— 2j} 

Un,2n-i=«0 2^ ' ( 31 ) 

A:— max{0;n— 2j'} 
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In particular, the leading asymptotic term of these coefficients as n — > oo and j 
fixed is given by 

jn-2j fj/\2(n-2j) j 
Un.an-J = «0Pi^ g»- 2j - (n _ m (1 + O(-)), 

(•_ •- ) n-2j/ n2(n-2j) 1 

t^an-j = 8 ;- 2i( ;: 2i) , (1 + O(-)). (32) 

Thus we have reduced the problem of the evaluation of the asymptotics of 
the coefficients u n ^ n -j v n,2n-j for large n to the computation of the leading 
coefficients of the polynomials Pj(z), Qj{z). In fact, using (|24| ) or ( p6|) and (p3|), 
it can be shown that the coefficients Pj,2j , <Zj,2j satisfy the recurrence relations 

(z/)V it/ 

Pj,2j = -l 2S . Pj-l,2(j-l) + 2~ 2^ Pk,2kPl,2iq m ,2m + 

k + l+ m =j 

H 2 6 ? 2 ^J-b2(j-l) ^ ?7-l,2y-l)J - 2^ Pk,2k(Pl,2l ~ qi,2l)Qm,2m, 

J J k,l,m=0 

k + l + m=j 

?j,2j - » 25 . gj-l,2(j-l) - 2^ 2^ Pk,2kqi,2iq m ,2m ~ (•«) 

^ ^ fe,l!,m=0 

k + l + m=j 

^2 W-1.2(j-l) - <7j-l,2(j-l) J + Jwj 2^ Pk,2k\Pl,2l - qi,2l)q m ,2m- 

fc + l + m=j 

Similarly, the coefficients it nj o, w «,o for the non-logarithmic terms appears 
from (|3l| ) for j = 2n, and are given simply by 

Unfi = U p 2 n,n, «n,0 — V Q q 2n ,n ■ (34) 

Thus the problem of evaluation of the asymptotics of the coefficients u n fi, Vn,o 
for n large is equivalent to computation of the asymptotics of the minor coeffi- 
cients in the polynomials Pj(z), qj(z). However, the last problem does not allow 
a straightforward solution because, according to (H), the sectional generating 
functions for the coefficients u n , v n are given by 

oo oo 

F(z, 0) = J2 u n,oz n , G(z, 0) = J2 W\ 

n=0 n=0 

and solve the separated Nonlinear Schrodinger equation 

+ izF z = vF + \zF" + F 2 G, 
8 

iG c -izG z = vG+lzG" + FG 2 . (35) 
8 



3 Discussion 



Our consideration based on the use of generating functions of different types 
reveals the asymptotic behavior of the coefficients u n 2n-ji v n,2n-j as n — > oo 
and j fixed for the long time asymptotic expansion (j^) of the generic solution 
of the sNLS equation (|l|). The leading order dependence of these coefficients on 
n is described by the ratio ( n ° 2 j)! • 

However, the behavior of the coefficients u Ut k, corresponding to other 
directions in (n, fc)-plane can be dramatically different. For instance, the gen- 
erating functions for u ni o, v n .o solving equation (]3^) are related to the so-called 
isomonodromy solutions [^o) of the system (0). These solutions are character- 
ized by the formal asymptotic expansion (E) without log-terms. Actually, it 
takes place when parameter v = —uqVq in (2) is independent of Ao = —x/2t. 
Concerning the asymptotic behavior of the coefficients u„.o, v n fl as n — ► oo, it 
is the plausible conjecture that tt n> o, w n ,o ~ c n T(^ + d). The investigation of the 
Riemann-Hilbert problem for the sNLS equation yielding this estimate will be 
published elsewhere. 
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